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USE OF ARBITRARY QUASI-ORTHOGONALS FOR CALCULATING FLOW DISTRIBUTION 
ON A BLADE-TO-BLADE SURFACE IN A TURBOMACHINE 


by Theodore Katsanis 
Lewis Research Center 

SUMMARY 

A method, suitable for computer programing, of analyzing flow on a blade -to -blade 
surface of a turbomachine is presented. The method is based on an equation for the ve- 
locity gradient along an arbitrary quasi- orthogonal between blades and is similar to a 
method using quasi-orthogonals in a meridional plane. With this method, a streamline 
analysis can be made for any blade-to-blade stream surface. This surface, if desired, 
may be assumed to be a surface of revolution generated by a meridional streamline ob- 
tained from a meridional streamline analysis. On this stream surface a two-dimensional 
solution for the velocity and pressure distributions is obtained. With several such blade- 
to-blade solutions, the velocity distribution throughout the rotor passage can be calcu- 
lated. Simplifying assumptions for upstream and downstream conditions are made for 
the purpose of readily obtaining a reasonable approximation near the inlet and outlet. 

The method was applied to a radial inflow turbine with splitter vanes. These results 
are given as a numerical example, and the blade surface velocities obtained are compared 
with those obtained by a simple linear velocity method. The comparison shows a large 
discrepancy between the linear velocity method and the method of this report near the 
outlet, but most of the rest of the blade is in reasonable agreement. The Fortran cpnw 
puter program is included. I 

r 

INTRODUCTION 

Quasi-three-dimensional methods have been developed for analyzing flow through 
mixed-flow turbomachines. One method is based on a meridional plane solution, with 
blade surface velocities calculated by some approximate formula. In reference 1 a 
method is given for obtaining a meridional plane solution using streamlines and quasi- 
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orthogonals to establish a grid for the solution. A quasi -orthogonal is defined to be any 
curve that intersects every streamline between the flow boundaries exactly once, as does 
an orthogonal to the streamlines. A quasi-orthogonal, however, is not necessarily or- 
thogonal to any streamline. The usefulness of quasi -orthogonals lies in the fact that they 
can be predetermined in some convenient manner, and they remain fixed regardless of 
any change of streamlines. Using this technique it was possible to develop a computer 
program that would calculate the velocity and pressure distributions without any interme- 
diate graphical procedures or hand calculations even for turbomachines with wide pas- 
sages and a change in direction from radial to axial within the rotor blade. 

From the meridional solution it is possible to obtain blade surface velocities by sev- 
eral methods. One method is by means of an approximate formula based on the assump- 
tion of linear velocity between blades and absolute irrotational flow (ref. 2) . This gives 
good results within the blades away from the inlet or outlet. Another method is the use of 
potential flow equations that may be solved by relaxation methods (ref. 3). The relaxa- 
tion method gives good solutions but is a rather tedious computation, even using a com- 
puter. It was realized that the basic idea used in reference 1 to obtain a meridional so- 
lution could be applied to obtain a blade -to -blade solution. By extending the solution up- 
stream and downstream of the rotor, a reasonable solution resulted throughout the rotor, 
and the solution was obtained with a reasonable effort. 

This report presents the analysis method for a blade -to -blade solution using quasi- 
orthogonals and discusses the numerical techniques required for obtaining solutions 
using a digital computer. The computer program developed for the blade -to -blade anal- 
ysis during this study is included. As a numerical example of the application of the anal- 
ysis method, a radial-inlet mixed-flow gas turbine with splitter blades and of high spe- 
cific speed is analyzed. 


METHOD OF ANALYSIS 

The blade -to -blade analysis is based on the assumption of a stream surface extend- 
ing from the pressure surface of a complete blade to the suction surface of the next com- 
plete blade, with the possibility of having a splitter blade between the complete blades. 
Although not necessary, it will be assumed that the stream surface is a surface of revo- 
lution. The shape of the surface of revolution may be determined by a hub-to-shroud 
streamline analysis, as indicated in reference 1. A separate blade -to -blade solution 
could be obtained for each streamline in the meridional solution. In many cases three 
blade-to-blade solutions would be adequate, one at the hub, one at the 50-percent meri- 
dional streamline, and one at the shroud. In the analysis presented herein, a correction 
for a loss in relative total pressure is included in the continuity equation to account for 
blade losses. 
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To obtain a reasonably good solution to the equations near the inlet or outlet of the 
rotor blade, it is necessary to extend the solution beyond the blade itself. The location 
of the streamline that passes through the stagnation point on the leading or trailing edge 
of the blade is not known. This location, however, can be approximated by various 
methods to obtain a better solution near the inlet or outlet. The basic method for a 
streamline analysis on a blade-to-blade surface for any type of turbomachine rotor or 
stator is given in the next section. This is followed by the Radial Flow Analysis which is 
applicable only to radial turbines. 


General Flow Distribution Analysis 


Analytical equations . - An equation for the directional derivative of the relative ve- 
locity based on an assumption of nonviscous flow is derived in reference 1. This differ- 
ential equation is used to calculate the flow distribution by numerically integrating the 
equation along an assumed quasi-orthogonal. The directional derivative of the velocity 
as given by equations (B13) and (B14) of reference 1 is 


dW dr .dz de i f 6 * 1 i 
— = a — + b — + c — + — I — 

dq dq dq dq W \ dq 



( 1 ) 


where 


a - W ros2 g cos a + sta a cos 0 -3L . 2a, si, p 

r_ r dm 


2 dW 

, W cos /3 sin a „„ „ a m 
b — + cos a cos p 

r dm 


/dW 

c = W sin a sin (3 cos /3 + r cos j3( + 2a> sin a 

dm / 


( 2 ) 


The coordinate system and notation are shown in figure 1. (All symbols are defined in 
appendix A. ) 

These equations give the velocity gradient along an arbitrary quasi-orthogonal on a 
stream surface where q is the distance along the quasi-orthogonal and where W is 
considered as a function of q alone. Note that q may be considered as a function of d, 
hence W may be considered a function of 6 alone, with the understanding that W is to 
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be on the quasi- orthogonal. Using the fact 
that 



Figure 1. - Coordinate system and velocity components. 


dW _ dW dq 
dd dq de 

results in the following relation being ob- 
tained from equation (1) 


dW 

de 


dr , dz 
a— +b— + c 


de 


de 


(3) 


along the quasi- orthogonal since 
dhj/de = dx/de = 0. Note that equation (3) 
is valid even if the stream surface is not a 
surface of revolution. 

In the examples given in this report, 
the blade -to- blade stream surface is as- 
sumed to be a surface of revolution, so 
that the quasi-orthogonal can be chosen a 
circular arc as shown in figure 2 with 



Figure Z - Quasi-orthogonal on surface of revolution. 


dr 

de 


^=0 

de 


An initial value must be determined 
for the integration of equation (3). This 
value is determined by the fact that the 
weight flow between blades is known. The 
weight flow across a fixed line from blade 
to blade can be computed by 


w = N 



An r d0 


(4) 


where An is the normal distance between 
two closely spaced streamlines in the 
meridional plane and can be determined 


from a meridional plane streamline analysis, as indicated in reference 1. The weight 
flow w is that portion of the total weight flow which is calculated to flow between the 
streamlines that are used to determine An. The density is calculated by the following 
equations, which are from equations (3) to (5) in reference 1 


P = 



y/(r- 1) 




where 


and 



W 2 + 2o)A - co 2 r 2 
2c p t ; 


( 5 ) 


( 6 ) 



T! 


2a)A - co 2 r 2 
2c T! 

pi 


( 7 ) 


Equations (5) to (7) and W m = W cos /3 give the numerical data for equation (4), which 
can then be numerically integrated. 

Numerical techniques and procedure . - The parameters a, f3, r , dW m /dm, and 
dW^/dm must be evaluated for use in equations (2) and (3). For this, a streamline 
geometry must be established. First, straight lines from blade to blade are established. 
These lines are the quasi- orthogonals. Since the assumed stream surface is a surface 
of revolution, the quasi- orthogonals can be chosen so that r and z are constant. For 
an initial approximation, each quasi-orthogonal is divided into a number of equal spaces. 
The parameters can now be evaluated in essentially the same manner as in reference 1, 
using spline fit curves for approximating first and second derivatives (see appendix C of 
ref. 1). The angle P is given by 


tan |3 = r — (8) 

dm 

where d0/dm can be calculated from a given streamline pattern by the use of a spline 
fit curve. 

After the parameters in equation (3) have been evaluated, the equation can be inte- 
grated by a Runge-Kutta method such as indicated by equation (13) of reference 1. Con- 
tinuity is used to determine the initial value of W in the numerical integration of equa- 
tion (3). After equation (3) has been integrated, the resulting flow distribution can be 
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used to obtain an improved streamline pattern. The details of this procedure is given in 
reference 1 in the section Numerical Techniques and Procedures. 


Radial Turbine Analysis 

As an illustration of the use of the general flow distribution analysis, the method will 
be applied to a radial inflow gas turbine with splitter blades. The specific considerations 
for the analysis of this type of turbine are given herein. 

Upstream analysis of radial turbine rotor . - In this section, equations are developed 
that are useful in extending the solution upstream for a rotor with radial inflow. At the 
inlet to the rotor of a radial inlet turbine, the relative flow direction usually is not paral- 
lel to the mean blade surface but is at a negative angle to the blade (i. e. , Wg is opposite 
to the direction of rotation). This means that inlet conditions cannot be uniform from 
blade to blade, since the flow direction must be parallel to the blade at the blade surface. 
For the purposes of the analysis, the solution was extended upstream one station where 
uniform inlet conditions were assumed. In most cases inlet conditions are specified at 
the inlet to the rotor. These conditions will be taken to be average from blade to blade. 
From this, the conditions at the upstream station can be calculated. 

■ From the assumption of free vortex flow, d(rV^)/dt = 0. This means that 
rV a = r.V fl • = A where the subscript i is used to designate a reference point such as 

U 1 ( 7 , 1 

the tip of the rotor. Also, by continuity pr W m An = p.mW m . An., so that 

W 

m, i _ rp An 
W m Vi An i 


These relations and the fact that tan /3 = W^/W m are used to calculate /3 as follows 


tan /3 = tan j3. 


W™ i W, 

m, l 9 


W M W m 


= tan 0. — ~ Vg 


lwr i- V e,i r iPi An i 


This gives 


tan /3 = tan /3. — - P An . 

1 corf - a Pi An i 


( 9 ) 
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For inward flow dr/dm = -1, since 
dm is positive in the direction of flow 
(inward) , but dr is positive in the out- 
ward direction. Since tan /3 = r d0/dm, 
this gives d0/dr = -(tan 0)/r. 

For the equation of the relative 
streamline, integrate 

de _ tang PAn tang, A \ 

dr r 'p, - *)V / 


Figure 3. - Relation between stagnation and free vortex relative streamline. 

This expression can be integrated by 

n . 

parts if u = p An and v = X log r - (wr^/2) in the formula /« dv = uv - f v du. This 
gives 


e - 9; 


tan /3. 


i,F 


( wr i - in i 


p An IX log r 



Pi An i| 


log 


r i 




for free vortex flow. 

If the last integral is approximated by the trapezoid rule, the following equation is 
obtained 


9-9 


h 



(p An + p. An^ 
2p. An t 



From figure 3 it can be seen that 9 - 9^ can be expected to be slightly less than 
9-9 ■ w , therefore it is assumed that 9 - 9- - M(0 - 0. F ), which gives 

i ? r 1 l, r 


9 - 0. = M 



(p An + p i Anj) 
2 Pi A^ 



( 10 ) 


where M is expected to be slightly less than 1. 
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For purposes of computation, a value of M is assumed, and after a blade- to-blade 

solution has been computed, the average value of V„ . between the blades is checked. 

y , 1 

If this average is not sufficiently close to the specified value of V fl ., a different value of 
M can be tried. It was found that the precise value of M was not critical in its effect 
on the solution. 

To completely specify a spline fit curve passing through a set of points, it is neces- 
sary to specify a condition at each end of the curve. In using the spline fit for approxi- 
mating a streamline in reference 1, it was assumed that the streamline was straightening 
out at the end, that is, the second derivative at the end point was specified to be one-half 
the second derivative at the next point. For the blade-to-blade solution at the inlet, how- 
ever, this is not a reasonable assumption. The coordinates used to specify the stream- 
line are 9 and m. For a radial inlet, d0/dm = (tan j3)/r, where tan 0 is given by 
equation (9). Thus, the derivative can be specified at the first point. This gives a better 
approximation to the streamlines. 

A similar situation occurs in calculating dW^/dm by means of a spline fit curve. 

In this case, 


dW, 

dm 


dV fl r i V fl i 

— 2 +aJ = 1 

dr 2 


+ o> 


for free vortex flow. Hence, at the first upstream station, 


dW, 

dm 



( 11 ) 


where r is the radius at the upstream station. 

For the streamlines that lead to the blade surfaces, there is a sharp curvature just 
before the blade followed by almost no curvature on the blade itself. To specify this 
streamline more accurately, additional points on the blade surface were specified be- 
tween the leading edge and the next station. Uniform inlet conditions are assumed at the 
upstream station. For the downstream end of the streamline it is still reasonable to as- 
sume the streamline is becoming straight. 

Downstream analysis for rotor with axial outlet . - It is assumed that the streamline 
downstream from the stagnation point on a rotor with axial flow at the outlet has a con- 
stant angle 0. This downstream angle of 0 can be determined from the blade angle at 
the outlet by correcting this angle for blade blockage. In the solution of the problem, this 
assumed streamline downstream from the blade is treated as if it were an extension of 
the blade. If the assumed streamline is close to the actual streamline, the velocity will 
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Figure 4 - Hub-shroud profile with streamlines used for blade- to- blade analysis. 


have little variation downstream of the rotor. If the variation is too great, it can be re- 
duced by adjustment of the assumed streamline, however, it has been found that a small 
variation in the assumed streamline has very little effect on the solution within the rotor. 

Analysis of rotor with splitter blades . - The method of analysis presented herein 
has been applied to a turbine rotor with splitter blades. It can be reasonably assumed 
that the weight flows through the passages on either side of the splitter blade are equal. 
The 50 percent streamline from blade to blade, therefore, meets the stagnation point on 
the splitter blade. A detailed analysis of the streamline shape in the neighborhood of the 
leading and trailing edges would be difficult and has not been attempted. Instead, the 
streamlines representing the suction and pressure surfaces of a blade are assumed to 
join one station beyond the end of the blade. This does not give detailed information on 
what is happening in the immediate neighborhood of the leading or trailing edge, but has 
little effect on the solution a short distance away. In the solution of the equations, the 
weight flow on either side of the blades is assumed to be 50 percent of the total weight 
between complete blades. This is used to evaluate the initial value for the integration of 
equation (3) on either side of the splitter blade independently. Downstream of the split- 
ter, equation (3) is integrated across the entire width of the passage. 


NUMERICAL EXAMPLE 

The method outlined has been applied to the analysis of a small radial flow turbine. 
The calculations were made on a digital computer. The hub- shroud profile is shown in 
figure 4. The blade has radial elements, except near the trailing edge of the splitter 
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Relative angle 


V 



Distance along meridional streamline, m, in. 
(c) Shroud. 

Figure 5. - Blade-to-blade streamline pattern. 
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blade, where the taper is not the same on both sides of the blade. The pertinent data for 
the case analyzed is given in the following list: 


Total number of blades 22 

Number of splitter blades 11 

Inlet radius, r^, in 3. 01 

Rotational speed, rpm 38 500 

Fluid Argon 

Weight flow, lb/sec 0.611 

Inlet total temperature, T|, °R 1950 

Absolute tangential velocity at inlet, V 0 ^ , ft/ sec 883 

Inlet relative flow angle, /3., deg -25. 1 

Inlet total pressure, p!, lb/sq in. abs 12. 97 

Loss of relative total pressure between inlet and outlet (assumed to vary 

linearly along streamline), Ap", lb/sq in 0.25 


The normal blade thickness was given by means of tabulated values on a grid. Blade 
thickness at any given point was obtained by linear interpolation. It was assumed that 
h! and A were constant across the inlet. 

The first step in the solution was to obtain a meridional plane solution as described 
in reference 1. This solution was extended downstream by assuming a mean stream sur- 
face downstream of the blades. Twenty meridional streamlines were calculated. The 
meridional streamline pattern obtained at the hub, mean, and shroud is shown in fig- 
ure 4. This gives the coordinates for a blade-to-blade surface at the mean meridional 
streamline, as well as the streamline spacing on the quasi- orthogonal s at hub, mean, 
and shroud. With this, the normal streamline spacing and the blade coordinates can be 
calculated at each of the three surfaces. The blades are extended by streamlines up- 
stream and downstream as explained in the sections on upstream and downstream analy- 
sis of a rotor. A computer program has been written for a combined meridional and 
blade-to-blade analysis, in which all the input data for the blade-to-blade analysis are 
computed by the program. The combined program is available to qualified persons if a 
written request is addressed to the author at the Lewis Research Center. The part of 
the Fortran computer program used for the blade-to-blade analysis for the numerical ex- 
ample is given in appendix B. The solution obtained by the computer program gives 
blade-to-blade streamlines in 6 - and m-coordinates. These streamlines have been 
plotted in figure 5. 

Figure 6 shows the blade loading on main and splitter blades at hub, mean, and 
shroud. The velocities are increasing over most of the blade; however, there is indi- 
cated a large negative velocity gradient on the pressure surface at the inlet near the hub. 


11 



Absolute velocity, V, ft/sec Relative velocit ''' W ' ft,sec 





(c) Shroud. 


Figure 6. - Blade velocity distribution for numerical example. 



Fraction of distance along meridional streamline 


Figure 7. - Shroud absolute velocity distribution for numerical 
example. 


It is even more severe at the end of the 
splitter blade near the shroud. This 
indicates the splitter should be ex- 
tended further from the aerodynamic 
point of view; however, the splitter is 
already longer than needed at the hub, 
and stress considerations may not per- 
mit making the splitter extend beyond 
the point where it ends on the hub. 
Figure 7 shows the absolute velocity 
distribution on the stationary shroud. 
Here, there are large negative velocity 
gradients. Of course, a negative ve- 
locity gradient cannot be eliminated 
entirely with the outlet absolute veloc- 
ity lower than the inlet absolute veloc- 
ity although it can be minimized by 
careful design. On the pressure sur- 
face, the nearly linear distribution is 
about the best possible for avoiding 
flow separation, under these circum- 
stances. However, along the shroud 
there is an increase, and then a 
greater decrease in a shorter dis- 
tance. It is difficult to avoid some- 
thing like this if there is to be any 
loading of the blades at the shroud. 

This type of analysis could be very 
useful as a design tool, since it points 
up the location of possible flow separa- 
tion. Modifications can be made in 
the geometry to improve on the veloc- 
ity distribution, until a good design is 
evolved. 

It is interesting to compare the 
blade surface velocities obtained here 
with those obtained by the method of 
reference 2, which is based on the as- 
sumption of linear velocity variation 
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(c) Shroud. 


Figure 8. - Comparison of blade surface velocities with those obtained 
by linear velocity method of reference 2. 

between blades. The comparison is shown in figure 8. The linear velocity method gen- 
erally shows slightly greater blade loading over most of the blade. Near the outlet and 
just past the splitter at the shroud there is a rather large discrepancy between the linear 
velocity method and the blade-to-blade streamline analysis. 


CONCLUDING REMARKS 

A method for determining the flow distribution in a blade-to-blade surface of a turbo- 
machine is presented. This method is based on a streamline analysis and can be used in 
conjunction with a meridional streamline analysis to obtain a quasi-three-dimensional 
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solution. The method was applied to a radial inflow turbine with splitter blades. The 
results are given as a numerical example, and the Fortran computer program is in- 
cluded. The blade- surface velocities obtained are compared with those obtained by the 
linear velocity method of reference 2. The comparison shows a large discrepancy be- 
tween the linear velocity method and the method of this report near the outlet but that 
most of the rest of the blade is in reasonable agreement. 

The basic ideas of the methods of streamline analysis on both blade-to-blade sur- 
faces and meridional planes are applicable to any type of turbomachine stator or rotor 
for either compressible or incompressible flow. The method appears to give reasonable 
results for the velocity distribution and for pressures throughout the passage. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 17, 1965. 
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APPENDIX A 




SYMBOLS 

a parameter, eq. (2) w weight flow between adjacent blade 

b parameter, eq. (2) to-blade stream surfaces 

c parameter, eq. (2) z axial coordinate 


c specific heat at constant pressure, 

(ft) (lb) /(slug) (°R) 

g acceleration due to gravity, 

ft/sec^ 

h enthalpy, (ft) (lb) /slug 

M coefficient, eq. (10) 

m distance along meridional stream- 

line, ft 

N number of complete blades 

An distance between adjacent stream- 

lines in meridional plane, ft 

p pressure, lb/sq ft 

Ap" loss in relative total pressure 

between inlet and any point, 
lb/sq ft 

q distance along arbitrary three- 

dimensional curve, ft 

R gas constant, (ft) (lb) /(slug) (°R) 

r radius from axis of rotation, ft 

r radius of curvature of meridional 

c 

streamline, ft 
T temperature, °R 

V absolute fluid velocity, ft/sec 

W relative fluid velocity, ft/sec 


a angle between meridional stream- 

line and z-axis, rad 

/3 angle between relative velocity 

vector and meridional plane, 
rad 

y ratio of specific heats 

9 relative angular coordinate, rad 

A prerotation; r.V„ ., sq ft/sec 

1 C7, 1 

p mass density, slugs/cu ft 

& rotational speed, rad/sec 

Subscripts: 

F free vortex flow 

i inlet 

m component in direction of 

meridional streamline, ft 

p pressure surface of blade 

r radial component 

s suction surface of blade 

9 tangential component 

Superscripts: 

' absolute stagnation condition 

** relative stagnation condition 
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FORTRAN PROGRAMS 
Description of Main Program 

The Fortran program listed here is for the blade-to-blade analysis of a radial inlet 
turbine with the same number of splitter blades as complete blades. This program was 
used to obtain data for the numerical example. It is written in Fortran IV and was run 
on the IBM 7094 digital computer. A meridional plane solution is required as input for 
this program. Reference 1 contains a program for the meridional plane analysis. The 
blade-to-blade program listed here follows the procedure given in the section Numerical 
techniques and procedure . 

The list of program variables given later in this section indicates the equation which 
is used to calculate a variable or the one in which it is used. In the program, the number 
of the streamline is indicated by K and the number of the quasi- orthogonal by I. Num- 
ber 1 is the station upstream of the inlet or the suction blade surfaces. The program is 
written so that all linear measurements are in inches, angles are in degrees, and pres- 
sure is in pounds per square inch for both input and output. Units are changed to feet 
and radians for computation in the program. All other quantities are in the units speci- 
fied in appendix A. 

It will be noted that a complete listing of input data cards is printed out. In the sam- 
ple program, for example, the listing gives all the data used as input for the program. 

All input cards precede the comment card END OF INPUT STATEMENTS. 

Program Variables and Definitions 

A temporary storage 

AB(J) temporary storage 

AC(J) temporary storage 

AD(J) temporary storage 

AE(J) temporary storage 

A.L(I) a 

ALM 
AR 
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A (input variable) 
R (input variable) 



* 


B 

BA(K) 

BCDP 

BETA(I, K) 
BE TIN 
C 

CAL(I) 
CBETA(I, K) 
Cl 

CORFAC 

CP 

CURV(I, K) 
DELTA 
DENSTY 
DN(I) 

DTDM(I, K) 

DTDMI 

DWDM (I, K) 

DWDM1 

DWDT 

ERROR 

ERROR 1 

EXPON 

GAM 

HT 

I 

IM2 

IND 


temporary storage 

ii. 

total weight flow between suction surface of blade and K streamline 

integer (input variable) , 1 will give THETA and WA as output on cards in 
binary form after final iteration; 0 will cause this to be omitted 

0 

/ 3 . (input variable) 
temporary storage 
cos a 
cos /3 

stagnation speed of sound at inlet 

percentage of calculated streamline correction to be used for next itera- 
tion 

C P 

curvature of streamlines in m, 6 plane 
calculated streamline correction 
PS 

An (input variable), eq. (4) 
de/dm, eq. (8) 

de/dm at upstream station (input variable) 
dW 0 /dm, eq. (2) 

dWg/dm at upstream station, eq. (11) 

<3W/d0, eq. (3) 

maximu m calculated streamline calculation for present iteration 
ERROR from previous iteration 
l/(y - 1), eq. (5) 

Y 

a e 

subscript to indicate number of quasi- orthogonal, 1 is upstream of inlet 
1-2 

code number used by subroutine CONTIN 
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IND 2 

IP2 

ITER 

ITNO 

J 

K 

KA 

KHMX 

KHP1 

KHP2 

KMX 

KMXM1 

LAST 

MX 

MXBL 

MXP2 

MXSP 

MXSP1 

NPRT 

NULL 

PLOSS 

PRINT 

PRS(K) 

R(I) 

RC(I) 

RHO 

ROOT 

RUNO 


% 


code number used by subroutine CONTAL 
I + 2 

number of iterations to be performed after ERROR is less than TOLER or 
with an increase in ERROR 

iteration number 

subscript 

subscript to indicate number of streamline, 1 at suction surface and KMX 
at pressure surface 

subscript 

KMX/2 

(KMX/2) + 1 

(KMX/2) + 2 

number of streamlines 

KMX-1 

KHMX-1 

number of quasi- orthogonals 

number of quasi- orthogonal at end of complete blade (input variable) 

MX + 2 

number of quasi- orthogonal at end of splitter blade (input variable) 

MXSP + 1 

data are listed for every (NPRT) streamline (input variable) 
dummy variable, not used 
Ap M at outlet, eq. (5) 

logical variable, when PRINT = TRUE output is listed 
P 

r (input variable) 
temporary storage 
p! g (input variable) 

V~2 

integer, run number 
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RR 


integer, 1 indicates velocities and streamlines from the previous run will 
be used as the initial estimate for the present case and 0 indicates that 
uniform velocity and equally spaced streamlines will be used as the ini- 
tial estimate for the present case 

SA(I, K) sin a sin 0 cos / 3 , eq. (2) 

SAL(I, K) sin a 

SB (I, K) r cos /3 ^dWg /dm) + 2u> sin a] , eq. (2) 

SBETA(I, K) sin 0 

SM(I) distance from inlet along meridional streamline 

SMAL(I) distance from inlet along meridional streamline with additional points just 

inside rotor inlet to make spline fit curve approximate blade more 

closely 

SRW integer (input variable) that will cause subroutines to write out data for 

certain values, used in debugging; 0 indicates that none of the sub- 
routines will write out extra data 

T-CURV curvature of streamlines in 9- and m-coordinates 

TEMP T! (input variable) 

THAL(I, J) 6 - coordinates of blades with additional points just inside rotor inlet 
THETA(I, K) 9- coordinates of blades and bounding streamlines (input variable) 

TIP T/T} 

TOLER (input variable); if maximum calculated streamline correction is less than 

TOLER, iterations are considered to have converged and desired output 
is printed 

TPP1P T"/T! 

TYPE integer (input variable), 1 indicates THETA and WA are given as input on 

binary cards and 0 indicates these quantities will be calculated by pro- 
gram 

V V 

VTHETA V Q 

W to (input variable) 

WA(I, K) W 

WAS W* (see eq. (13) in ref. 1) 
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WASS 

W** (see eq. (13) in ref. 1) 

WM 

W m 

WT 

w (input variable) 

WTFL(K) 

calculated weight flow between assumed blade-to-blade stream sur- 

1L 

faces and between suction blade surface and K in streamline, eq. | 

WTH 

w/2 

WTHETA(I, K) 


WTOLER 

allowable tolerance on weight flow in satisfying continuity between 
blades (input variable) 

XN 

N (input variable) 

Z(I) 

z (input variable) 


Fortran Program Listing 


DIMENSION Z(21) ,R(21> ,DN(21 ) , SMI 21 ) ,BA<22) ,AB(22) ,AC(22) ,AL(21 ) , 
1RC121) , CAL (21 ) , SAL (21 ) ,PRS(22> , WTFL(22) , WTHE TA ( 2 1 , 22 ) ,DwOM(21 ,22) , 
2THETA( 21, 22 > »WA( 21,22) ,BETA( 21,22 ) ,SflETA( 21 ,22) ,CBETA( 21,22) , 

3SA ( 21,22) ,SB( 21,22) , CUR V< 2 1 , 22 ) , DTDM< 2 1 , 22 ) 

DIMENSION Zl(21 ,21 ) ,R1( 21,21 ) , VHH1 ( 21 ) , THHKHt 21) , THHKPI21 ) , 

1 T HHKMX( 2 1 ) , THM1 ( 2 1 ) , THMKH ( 2 1 ) , T HMKP (21) , THMKMXl 2 1 ) , T Hi, 1(21) , 

2 I HSKHt 21 ) ,TnSKP(21>, THSKMX121) 

01 MENS ION THAL( 23,4) ,SMAL( 23 ) , ALH 22 ) , At ( 22 ) ,UtNSTY( 22) ,bETA1{ 3) 
INTEGER RUNO, TYPE,U LOP, SR W,RR, HUB, SHROUD 
LOGICAL PRINT 
RUNU=0 

5 READ 15,1010)MX,KMX,MXSP,RR,W,WT,XN,GAM,AR 
I TNO 3 l 
KUNO=RUNO+1 
WRITE (6,1020) RUNG 

WRITE 16, 1010)MX,KMX,MXSP,RR,W,WT,XN,GAM,AR 
RE AO (6,1010) TYPE, 3C0P, SRW.MXRL, TEMP, ALM,RHO, TOLER, PLOSS 
WR 1TE(6,1010)TYPE, liCDP , SRW , MX&L ,TEMP,ALM, RHO »TOLER»PLOSS 
PL0S:> = PL0l>S*L44. 

READ (6,1010) NULL, NPRT, ITER, NULL, BET I N , WTOLER , CORFAC 

WRITE It., 1010) NULL, NPRT, ITER , NULL , BET I N, WTOLER , CDRFaC 

BETIn * BET IN/57. 20577 

KHMX - KMX/2 

KHP1 =• KHMX+1 

KHP2 =« KHMX+2 

READ (6,1030) ( THE TaI I , 1 ) , I = 1 , MX ) 

WR I T£ 1 1> , 1 030 ) (THETAf 1,1), 1 = 1, MX) 

RE AO (6,1030) ( THETA ( I , KMX ) , 1=1 ,MX) 

WRITE! fa, 1030) ( THETA ( I ,KMX) , 1=1 ,MX) 

READ (6,1030) ( THE TAI I , KHMX ) , 1*1 , MXSP ) 
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non 


WRITEib, 1030) f T H E T A ( I.KHMX), 1=1, MXSP) 

READ ( 5, 1030) ( THETA ( 1 , KHP1 ) , I =1 , MXSP ) 

WR l Tt l 6 » 1030 ) ( THETA( I, KHP1 ) , 1 = 1, MXSP) 

lF(RR.tQ.i) GO TO 90 

READ (5,1030) ( Z ( I ) , 1=1, MX) 

WRITE (6,1030) ( Z( I ) , 1=1, MX) 

READ ( 5, 1G30) I R ( I ), 1=1, MX) 

WR1 TE( t»» 1030) ( R ( I) , 1 = 1, MX) 

READ I 5,1030) (0N( I ) , 1 = 1, MX) 

WRITE! 6,1030) (DN( I) , 1=1, MX) 

00 9 1=1, MX 

Z(I)=Z( D/12. 

R( I )=R( I)/ 12. 

9 ON ( I ) =UN( 1 )/12. 

21 IFITYPt.EG.l ) GO TO 30 

WA ( 1, i)=WT/RHO/ON( 1)/RI 1 »/XN/( THETA ( 1, KMX )-THETA( 1,1) ) 

00 23 1=1, MX 

IF1 ( l.fcO.l ).UR.( l.GT.MXSP) ) THETA! I.KHMX1 = ( THETA( I ,KMXH-THETA 
1 (I,l))/2. 

IF((I. EO. 1 ). OR. (l.GT.MXSP)) THETA ( 1 ,KHP1 ) = THE TA ( I , KHMX ) 

DO 22 K=1,KHMX 

THE TA1 1 , K ) = FLOAT (K-l) /FLOAT! KHMX-DM THETA ( I.KHMX) -THETA! I, 1) ) 

1 ♦ THfcTAl 1,1) 

22 MA ( 1 , K) =MA (1,1) 

OU 23 K=KI(P 1 , KMX 

THtTA( I ,K) = FLOAT ( K— KHP 1 ) /FLOAT (KMX-KHP1)*(THETA( 1, KMX) — 

1 THETA! 1 , KHP1 ) ) +THETA! I ,KHP1 ) 

23 MA( I,K) = MA( 1,1) 

READ (5,1010) NEXT 

24 IF( NEXT .EQ.O) GO TO 45 

READ (5,1021) 1 1 ,K 1 , TH£TA( 1 1 , K 1 ) , 1 2 , K2 , THETA ( I2»K2 ) , 13 ,K3, THETA( 13, 
1 K3 ) , 14, K4, THETA( 14, K4) , 15, K5, THETA ( I 5, K5) .NEXT 

WRITE! 0,1021) 1 1 , K 1 , THETA ( 1 1 , K 1 ) , 1 2 , K2 , THETA ( 12 , K2 ) , 13,K3,THETA( 13, 
1 K3) « I4,K4,THETA( 14, K4), 15,K5,THETA( I 5, K5), NEXT 

GO TO 24 

END OF IMPUI STATEMENTS 

30 GALL BUREAU ( THETA! 1 , 1 ), THETA ( 21, 21 ) ) 

CALL tiCREAO I WA (1,1), WA(2l,21)) 

WRITE (6,1040) 

45 CONTINUE 

CP=AR*GAM/(GAM-i. ) 

Cl = SwRT ( GAM*AR*TEMP ) 

WRITE (6,1050) Cl 
KMXM1 = KMX-i 
MXP2 = MX-*- 2 
MXSP 1 = MXSP+1 
WTH = WT/2. 

EXPON = l*/( GAM-l *) 

TANtJTU = SIN(BETIN)/COS(BETIN) 

TANHT = TAN6IO«(ALM-W*R(l)**2)/(ALM-W*R(2)**2)*ON( 1)/DN(2) 

DTDM1 = T ANBT /R ( 1 ) 

DrtDMl = w+ALM/R ( 1 )**2 
R001 =S*RT ( 2 . ) 

SMI 1 ) = —SORT ( ( Z( 2 ) — Z ( 1 ) )**2+(R(2)-R( 1 ) )*»2) 

DO t.0 K = 1 , KMXM 1 

60 6 A (K)=ELOAT( K— 1 ) *WT/ FLOAT ( KMXM 1-1 ) 
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CALCULATE ALPHA AND SM 

DO 70 i = 1 » MX 
AtSI 1 ) =1 Zt I )-Kl I ) ) /KUOT 
70 ACI I) = l ZIIJ+RI I > ) / RUO T 

CALL SPL I NE ( AB, AC, MX, AL,RC) 

DO 80 1=1, MX 

AL ( I ) = ATAN ( AL ( I ) >-.78 5398 
CAL ( I ) = COS ( AL ( I) ) 

80 SAL (• 1 ) =S IN ( AL ( l > > 

DO 85 1 =2 ,MX 
J = 1 - 1 

8 5 SMI l)*SMI J )+S0RTt ( Z ( I ) -Z ( J ) ) **2 + ( R < 1 1 — R ( J) )**2) 

SMAU3) * SMI2)+.1*(SMI3)-SM(2)) 

SMALI4) = 5M(2)+.5*(SMI3)-SM<2) ) 

90 ERROR = 1000. 

CALCULATE BETA ON BLADE SURFACES 
BEGINNING OF LOuP FOR ITERATIONS 

91 DO 97 0=1,4 
K = 1 

lFlJ.fcw.2) K=KHMX 
I F ( J . ED . 3 ) K = KHP 1 
IFIJ.Ew.4) K = KMX 
DO 93 1=1,2 

THAU I, J) = THETA! I,K) 

93 SMAL ( I > = SMI I ) 

DO 94 1=3,4 

94 THAU I, J) = 1 HET A I 2 * K ) + 1 SMAL I I ) -SM I 2 ) ) / I SMI 3 ) -SM I 2 ) ) * I THET A ( 3 * K ) 
1 -THtTAI 2,K) ) 

DO 95 I =5,MXP2 
I M2 = 1-2 

THAU 1, J) = THETA! IM2.K) 

95 SMAL I I ) * SMI IM2 ) 

CALL SPL IN2 I >MAL , THAL I 1, J) , DTDM 1 , MXP2 , DTOMI 1 , K ) , AB > 

DO 96 1=1,2 

CURVII.K) = ABI I ) / 12 . / I l.+IDTDMI I » K ) / 12 • )**2)**1.5 
BETA! 1 , K ) = aTANIRI i)*DTDMl I, K) ) 

SBE1 A ( I , K ) = SINIBETA! I,K) ) 

C8E1 At i , K ) = C0S1 BETA! 1 , K) ) 

96 CONTINUE 

DU 97 1=3, MX 
IP2 = 1+2 

CUR V I I , K ) = aB{ IP2)/12./< l. + IDTDMI 1P2,K)/12. )**2>**1. 5 
bETAI1,KJ = ATANIRI I )*DTDMI IP2,K) > 

S8ETAI I ,K> = SINIBETAl I , K) ) 

CBETAU.K) = COSI BETA! I,K) ) 

97 CONTINUE 

PRINT = I I TEK.LE.O ) .OR. ( ITNO.LE.NC ) 

IF (PRINT) WR1TEI6, 1060) ITNO 

ERROR UERRCJR 
ERROR=u • 
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START CALCULATION UF FAkAMETERS 

LAST = KHMX— i 
UO 100 K=2,LAST 

CALL SPLIN2I SM, THETAl 1,K) ,0TUM1 ,MX,UTDM( l ,K) ,AB) 

DO 100 1*1, MX 

CUKV( 1 ,K) = aB< l)/12./Il.-MDTUM(I,K)/12.)**2)**i.5 
BE 1 A( 1 ,K)=ATaN(R( IMDTOMI l f K) ) 

SBETAl 1,K)=S1N( BETA I I,K) ) 

CBETAI I » K ) =CUS ( bETA( I,K) ) 

100 CONTINUE 

DU 101 K=KHP2 , KMXM 1 

CALL SPLIN2IsM,THETA( 1 , K ) , OTOfll , MX, IJTDM ( 1 ,K ) » Aft) 

DO 101 1=1, MX 

CURVI 1 ,K1 = aB(I)/ 12./( l.+IOTDMl I,K)/12.)**2)**1.5 
BETA! I,K)=ATaN(R( I)*UTUM{ I,K) ) 

SBLTAI 1,K)=S1N( BETA! I ,K) ) 

CBETAI 1 ,K ) = COS( BETA ( I , K ) ) 

101 CONTINUE 

UO 110 K* 1 » KMX 

00 1 0 5> 1 = 1, MX 

WTHETAI I , K ) =wA ( I,K)*SBET a( I,K) 

105 CONTINUE 

CALL SPLIN2I SM.WTHETAl 1 , K ) , OWOMl , MX , DWOM( 1,K) , AC) 
DO 110 1=1, Mx 

SAII.K) * SALI I )*SBETA( I,K)*C8ETA( I ,K) 

SB I I , K) * CBfcTAI I,K)*RI I ) * < 2 . *W*SAL < I )*OWDM( I,K)) 
110 CONTINUE 

END OF PARAMETER CALCULATION 

CALCULATE VELOCITY Dl STRIBUT ION, CHECK CONTINUITY 

DO 200 1*1, MX 

1 ND= 1 
INU2 = 1 
GO TO 130 

120 WAI I » 1 ) =• 5* WA ( 1,1) 

GU TU 130 

125 WAI I, 1 ) =2 . * WA ( 1,1) 

GO TU 130 

126 WA ( I , KHP1 ) = . 5* WA ( 1 , KHP 1 ) 

GO TU 130 

127 WAI I , KHP1 ) = 2.*WA( l.KHPl) 

130 IF! ( I . GT. I ) .AND. ( I.LE.MXSP) ) GO TU 1A2 
DO 1 AO K*2 , KMX 
J = K-l 

I F ( K. Eu.KHP 1 ) W A ( I , K ) = wA(l,J) 

IFIK.Ey.KHPl ) GO TU IAO 
HT = THETA( I, K) -THETA I I,J) 

WAS = wAI 1 , J ) + ( WA ( 1 , J ) *SA( I , J ) +S8I 1 , J ) )*HT 
WASS = WAI I, J )+I WAS*SA( I,K)+Srt( I,K) )*HT 
WAI I , K ) * ( WAS+WASS )/2. 

1 AO CONTINUE 
GU TO 1A8 

1 A2 DO 1AA K=? , KHMX 
J=K— 1 

HT=THfcIAI I, K) -THETA I I, J) 

WAS = wAl I,J)+( WAI I,J)*SA( I,J)+SB( I,J) ) *HT 
WASS = WAI I , J ) + ( WAS*SA( I ,K ) +SBI I,K) )*HT 
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WA4I,K) * ( WaS+WASS ) /2* 

144 CONT 1NUE 

l)G 1 4b K=KHPc,KMX 
J = K-i 

HT=THETA1 I ,K)~THETAI I, J) 

WAS = wA I I * J ) + ( WA ( I , J ) *SA I I , J ) +SB I I t J ) ) *HT 
WASS = WA< I,J)MWAS*SAl I,K)+SB( I,K) )*HT 
W A ( I , K ) * (WaS+WASSI/2. 

146 CONTINUE 
148 CONTINUE 

00 150 K*i,KMX 

TIP= i.-(WA( 1,K)**2+2.*W*ALM-(W*RI I) ) **2 ) /2./CP/TEMP 

IF( (-T1P.LT. .0) .AND. ( l.GT. 1) .AND. I I.LE.MXSP) . AND. ( K . GE . KHP 1 ) ) 

1 GO TO 126 
IFIT1P.LT. .0) GO TO 120 

TPP IP - 1.- 12.*W*ALM-(W*R( I ) > **2 ) /2. /CP/TEMP 

DENST Y I K ) = T lP**EXPUN*kriO-< T1P/TPP1P) **E XPON*PLOSS/AR/ TPP 1 P/ TEMP 

1 *32.1 7*SM( I ) /SMIMXBL) 

PRblK) = DENiTY ( K ) *AR*T lP*TEMP/32. 17/144. 

WM=WAt I , K ) *CbET A ( I , K ) 

A B 4 K ) = UENSTY ( K ) *WM*DN ( I ) *R ( I ) *XN 

150 ACIK)=*Th£TA( 1,K) 

CALL 1nTgRLIAC»A6»KHMX»WTFL) 

1 F I WTFL I KHMX ) • LE • • 0 )• GO TO 125 

CALL lNTGRLl aCIKHPI ) , ABI KHP1 ) , KHMX , WTFL 1 KHP 1 ) ) 

IFI I'WTFLIKMX) .LE. .0) .AND. I l.GT. 1) .ANU. I I.LE.MXSP) ) GO TO 127 

IFIWTFLI KMX ) . Lfc . . 0 ) GO TO 125 

IF{( I.GT.ll.ANO.I I.LE.MXSP) ) GO TO 155 

00 152 K*KHP 1 , KMXM 1 
AC(K) = AC ( K+ 1 ) 

152 WTFLIK) * WTPL(KHMX)+WTFL(K+1) 

IF I AB* I NT- WTFL CKMXM1) l.LE.WTOLER) GO TO 160 
CALL CUNT INI wA ( 1,1) , WTFL I KMXM 1 ) , I NO, I , WT ) 

1 F ( INO .N£ «6 ) GO TO 130 
GO TO 160 

155 IF (ABS(WTH-wTFLIKHMX) l.LE.WTOLER) IND=6 
IF LABS! WTH-wTFL ( KMX ) ) . LE . WTOLER ) IND2=6 
IFIIN0.NE.6) CALL CONT I N ( WA ( I , 1 ) , WTFL I KHMX ) , INO , I , WTH ) 

IFI IN02.NE.6) CALL CONT AL ( WA ( I , KHp 1 ) , WTFL l KMX ) , I N02 , I , WTH) 

I FI ( IN0.NE.6) .OR. ( IND2.NE.6M GO TO 130 
160 IFK I. GT.1 >. AND. I I.LE.MXSP) ) GO TO 165 
CALL SPLINT I WTFL , AC , KMXM1 , BA , KMXM1 , A3 ) 

GO TO 166 

165 CALL SPL1NTIWTFL, AC, KHMX, 0A, KHMX, AB) 

CALL SPL INTI WTFL I KHP 1 ) ,AC( KHP1 ), KHMX, BAl 1 ) , KHMX , AB I KHP 1 ) ) 

166 CONTINUE 

IFK I. GT.l). AND. I I.LE.MXSP) ) GO TO 168 
DO 167 K* 1 , KHMX 
KA = KmX-K+1 
J = KA- 1 

167 ABIKA) = ABIJ) 

168 00 170 K- 1 , KMX 
OELTA=»aBSIAB(K)-THETA( I,K) ) 

1 70 IFIDELTA.GT.tRROR)ERROR=DELTA 
IFI. NOT. PRINT) GO TU 178 
A = $Ml I ) * 12 . 

C = AL< 1 )*57. 29577 

D = Rl I ) * 1 2 . 

E = Z I 1 >*12. 
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on o 


F = UN( I)*12. 

WRITE (6,1080) I ,A,0,D,E,F 
WRITE (6,1070) 

00 176 K* 1 , KHMX , NPRT 
B=BETA( I,K)*b7.29577 
C = WTHET A ( I , K ) 4-W*R ( I ) 

WM = WA( I,K)*C8ETA( I,K> 

V = SWRT{ C**2+WM»*2) 

DWDT =* WA( I , K ) *SA ( I , K ) + Sti( I ,K ) 

WRITE (6,1090) THE T A( 1 ,K) , CURVI I » K ) ,B * WA( 1 , K ) , WTHETA( I , K ) , C , WM, V , 

1 PRS(K) »DENSTY(K) ,DTDM( I , K ) , DWDM1 I , K ) , SA( I ,K) ,SB( I,K) ,DWD1 

175 CONTINUE 
WRITE (6,1095) 

00 176 K=KHP1, KMX, NPRT 
B=BET A ( I , K ) *57 • 295 7 7 
C = WT HETA { l,K)*W*R( I) 

WM = mA( I»K)*CBETA( I,K) 

V = SURT( C**2+WM**2 ) 

DWDT = WA( I,K)*SA( I,K)+S3( I,K) 

WRITE (6,1090) THETA( I , K ) , CUR V ( 1 , K ) , B , WA { I,K) , WTHE TA ( I,K) ,C,WM,V, 

1 PRS 1 K ) ,DENSTY(K) , DTDMI I , K ) , DWDM ( I , K ) , SA< I ,K ) ,SB( I ,K) ,0MD1 

176 CONTINUE 
WRITE (6,1095) 

178 00 180 K*2 , KMXM1 

180 THETAt 1 , K > = ( 1.— CORF AC ) *THETA( I ,K > ♦CORFAC*AB ( K ) 

200 CONTINUE 

END OF VELOCITY CALCULATIONS. 

WRITE (6,1120) ITNO, ERROR 
I F ( ITER. LE. 01 GO TO 230 

I F ( (ERROR. GE.ERR0R1 ) .OR. (ERROR. LE. TOLER) ) I TER= I TER-1 
I T N0= 1 T N0+ 1 
GO TO 91 

230 IF(BCOP.NF.l) GO (0 5 

CALL BuDUMP (THETA(1,1),THETA(21,21) ) 

CALL BCDUMP ( WA(1,1), WA(2L,21)) 

GO TO 5 

1010 FORMAT (415, 6F l 0.4) 

1020 FORMAT (8H1RUN NO. 1 3, 10X, 25HINPUT DATA CARD LISTING ) 

1021 FORMAT (5(2I2,F8.5) , I 1) 

1023 FORMAT ( 32H THETA-CALCULATED AND/ UK INPUT ) 

1030 FORMAT ( 7F10.4) 

1040 FORMAT ( 10X24HBCD CAROS FOR THETA, WA ) 

1050 FORMAT ( 36HK STAG. SPEED OF SOUND AT INLET = ,F9.2,//> 

1060 FORMAT ( ///5Xl 3HI TEKAT ION NO. 13) 

1070 FORMAT ( 1 12 HK THETA T-CURV BETA WA WTHETA VTHETA 

1 WM V PRS DENSTY DTDM DWDM SA SB DW 

2DT ) 

1080 FORMAT ( 2X16H0UAS I-GR THUGONAL I 3 , 6X, 4HSM =,F7.4,9H ALPHA =,F7.2, 

1 5H R = , F 7 .4, 5H Z =,F7.4,6H UN =,F 7.4) 

1090 FORMAT ( IX , F9.4 , 2F8 .2, 5F8. 1 , F7. 2, F9. 5 , F7. 2 , F8. 0, F8 . 4, 2F9. 1 ) 

1095 FORMAT ( IH ) 

1120 FORMAT ( 18HJ ITERATION NO. 1 3, 10X, 24HMaX. STREAMLINE CHANGE = , 
IF 10.6 ) 

1200 F0RMAT1// 10X7HN0RMAL 14) 

1210 FORMAT. 1 7F 18 . t> ) 

END 
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Description of Subroutines 


The subroutines SPLINE, SPLINT, INTGRL, and CONTIN are described in refer- 
ence 1. The subroutine CONTAL is the same as CONTIN and is used so that a continuity 
calculation across a quasi- orthogonal can be computed on both sides of the splitter blade 
simultaneously. The subroutine SPLIN2 is based on the spline fit curve, with the deriva- 
tive specified at the first point of the curve. SPLIN2 gives the first and second deriva- 
tives at each point. The calling sequence for SPLIN2 is as follows: 

CALL SPLIN2(X, Y, YIP, N, SLOPE, EM) 


where 

X input array 

Y input array, function of X 

YIP, input, derivative at first point, dY/dX (X(l)) 

N input, number of X and Y values given 

SLOPE output array, first derivative, dY/dX 

2 2 

EM output array, second derivative, d Y/dX 


The subroutine SPLIN2 is as follows: 


SUbROUTINE SPLIN2(X,Y,yiP f N,$LOPE,EM) 

DIMENS ION X ( bO 1 , Y l 50) » S ( 50 ) » A ( 50 ) *B(50) » C ( 50 ) * F { 50) » W ( 50) » SB( 50 ) » 
10 C 50 > ,tM(50) » SLOPE ( 50) 

COMMON Q 
INTEGER Q 

00 10 1=2 tN 

10 S( I ) = Xl I >-Xt 1-1 ) 

N0=N- L 

00 20 1=2, NO 

a ( n=sm/ 6 . 

B< !>*(*( I )+SU + l> )/3. 

CID-SH + ll/6. 

20 E(n = (Y<l + n-Y(I))/G(U-i)-(Y(l)-YtI-l))/S(I) 

A ( N ) =— . 5 
till )=->! 2)/3. 

U(N)=i. 

C(l)=i(2)/6. 

F ( 1 ) = l Y ( 2 >-Y< 1) )/S( 2J-YIP 
F ( N ) =0 . 

Wll )=iU 1 ) 

^bi n=c( l ) / w 1 1 ) 
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G< i)=F{ 1)/W< i) 

00 30 i = 2 » N 

w( I )=B{ I)-A( 1 ) * SB ( 1-1 ) 

s«( n=u n/wi i) 

30 G( I ) = (FU)-A( I )*G( l-l) )/WC I) 

EM(N)=G(N) 

0U 40 1=2, N 
K=N+i— 1 

40 EM(K)=6<K)-S6<K)*EM(K+1) 

SLOPE l 1)=— S(2)/6.*(2.*EM( 1 )+EM(2) )+<Y(2)-Y(l))/S(2) 

0050 1=2, N 

50 SLOPE ( I) = i< I)/6.*(2.*EM( I ) +EM ( I — 1 ) ) + ( Y { I ) — Y ( I — 1 I )/S( I) 

IF (U.fcU,14> WRITE (6,100) N, ( X (I) , Y (I) , SLOPE (I) , EM(I) , 1= 1 ,N ) 

100 FORMAT (2X15HHO. OF POINTS =I3/10X5HX 15X5HY 15X5rtSL0PEI5X5H 
ltM /14F20.8)) 

RE TORN 
END 


Sample Output from Program 

The output given here is the listing for the case used in the numerical example for 
the mean blade-to-blade surface. It will be noted that there is an exact listing of all in- 
put data cards at the beginning of the listing. However, the value for WTOLER should be 
0. 00001 instead of 0. 0000 as rounded off in the listing. The input data listing is followed 
by the maximum calculated streamline change for each iteration, which is used as the 
criterion for convergence. After 26 iterations, there is convergence within the specified 
limit of 0. 001-radian maximum streamline change. At this time, streamline coordinates 
are printed out together with velocities, pressure, and certain parameters used in the 
calculations. The calculations are based on 20 stream tubes between the blades. Data 
are printed out only for every other streamline in this listing. 
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N: C O O C 


HUN NU. 1 iNPUT DATA CARD LISTING 

lb 22 7 G 4C3C.0000 C.0305 11.0000 1.6667 1245.0000 

G 0 12 1950.0000 221.5000 0.0247 0.0010 0.2500 

c 2 C -23.3600 C.0000 0.1500 

.4970 -0.5626 -0.5590 -0.5556 -0.5523 -0.5503 -0.5497 

.5552 -0.5626 -0.6674 -C.B503 -1.1139 -1.3724 -1.7284 

.0644 

0.0742 -C.0C86 -0.0122 -0.0156 -0.0190 -0.0207 -0.0222 

-0.0264 -0.0556 -0.1452 -0.3347 -0.5842 -0.8012 -1.1572 

-1.5132 

-0.2091 -0.^942 -0.2978 -C.3012 -0.3014 -0.3002 -0.3016 

-0.2091 -C.277C -0.2734 -0.2700 -0.2712 -0.2767 -0.2880 

0.3755 0.3555 G.4576 0.5551 0.7143 0.8720 1.0467 

1.2370 1.4521 1.7906 2.1000 2.3500 2.5500 2.8500 

3.1500 

3.3100 3.0100 2.6770 2.4028 2.1541 2.0036 1.8931 

1.8120 1.7460 1.6508 1.6262 1.5767 1.5606 1.5527 

1.5510 

0.0360 O.C378 0.0406 0.0429 0.0448 0.0457 0.0458 

0.0489 0.0522 0.0562 C.0602 0.0639 0.0658 0.0666 

0.0666 


STAG. SHfcfcO OF 

SOUND AT 

1NCET = 2011.55 




HERAT 1 ON 

NL. 

1 

MAX. STREAMLINE 

GHANGE 


0.047759 

1TEKAT1CN 

NO* 

2 

MAX. STREAMLINE 

CHANGE 

= 

0.066704 

ITERATION 

NO. 

3 

MAX. STREAMLINE 

CHANGE 

= 

0.056437 

ITERATION 

NO . 

4 

MAX. STREAMLINE 

CHANGE 

- 

0.046788 

ITERATION 

NO. 

5 

MAX. STREAMLINE 

CHANGE 

■= 

0.037940 

ITERATION 

NO. 

6 

MAX. STREAMLINE 

CHANGE 

= 

0.030659 

ITERATION 

NO . 

7 

MAX. STREAMLINE 

CHANGE 

X 

0.024828 

ITERATION 

NO. 

8 

MAX. STREAMLINE 

CHANGE 


0.020120 

ITERATION 

NO. 

9 

MAX. STREAMLINE 

CHANGE 

= 

0.016351 

ITERATION 

NO. 

10 

MAX. STREAMLINE 

CHANGE 

= 

0.013294 

ITERATION 

NO • 

It 

MAX. STREAMLINE 

CHANGE 

= 

0.010829 

ITERATION 

NU. 

12 

MAX. STREAMLINE 

CHANGE 

IS 

0.008846 

ITERATION 

NO. 

13 

MAX. STREAMLINE 

CHANGE 

= 

0.007244 

ITERATION 

NO. 

14 

MAX. STREAMLINE 

CHANGE 

= 

0.005955 

ITERATION 

NO. 

15 

MAX. STREAMLINE 

CHANGE 

S 

0.004932 

ITERATION 

NO. 

16 

MAX. STREAMLINE 

CHANGE 

= 

0.004065 

ITERATION 

No • 

17 

MAX. STREAMLINE 

CHANGE 


0.003352 

HERAT ION 

NO • 

16 

MAX. STREAMLINE 

CHANGE 

= 

0.002934 

11 ERA! ION 

Nu . 

15 

MAX. STREAMLINE 

CHANGE 

= 

0.002352 

I T EKAT I ON 

NO . 

20 

MAX. STREAMLINE 

CHANGE 

= 

0.002057 

Herat ion 

NO. 

21 

MAX. STREAMLINE 

CHANGE 

- 

0.001669 

ITERATION 

NU . 

22 

MAX. STREAMLINE 

change 

SS 

0.001507 

ITERATION 

NO . 

23 

MAX. STREAMLINE 

CHANGE 

= 

0.001218 

1 1 tRAT I ON 

NU . 

24 

MAX. STREAMLINE 

CHANGE 

= 

0.001095 

i T tKAT i ON 

No • 

2b 

MAX. STREAMLINE 

CHANGE 


0.000985 
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